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ABSTRACT 

Worked  performed  during  this  period  ineludes  the  investigation  into  eigenvalue  behavior  of  several  different  elasses  of  large  dimensional 
random  matriees.  They  are:  1)  a  elass  of  random  matriees  important  to  array  signal  proeessing 

and  wireless  eommunieations  with  the  goal  of  proving  exaet  separation  of  their  eigenvalues;  2)  an  ensemble  of  random  matriees  used  to 
estimate  the  powers  transmitted  by  multiple  signal  sourees  in  multi-antenna  fading  ehannels;  3)  another  ensemble  whose  eigenvalues  yield 
the  mutual  information  of  a  multiple  antenna  radio  ehannel,  for  whieh  a  eentral  limit  theorem  is  proven;  4)  ensembles  whieh  yield  robust 
estimation  of  a  population  eovarianee  matrix  with  applieation  to  array  signal  proeessing;  and  5)  a  sample  eovarianee  matrix  for  whieh  a  CLT 
is  studied  on  linear  statisties  of  its  eigenvalues,  whose  limiting  empirieal  distribution  of  its  eigenvalues  is  studied  with  applieation  toward 
eomputing  the  power  of  a  likelihood  ratio  test  for  determining  the  presenee  of  spike  eigenvalues  in  the  population  eovarianee  matrix. 
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Worked  performed  during  this  period  includes  the  investigation  into  eigenvalue  behav¬ 
ior  of  several  different  classes  of  large  dimensional  random  matrices.  They  are:  1)  a  class  of 
random  matrices  important  to  array  signal  processing  and  wireless  communications  with 
the  goal  of  proving  exact  separation  of  their  eigenvalues;  2)  an  ensemble  of  random  matri¬ 
ces  used  to  estimate  the  powers  transmitted  by  multiple  signal  sources  in  multi-antenna 
fading  channels;  3)  another  ensemble  whose  eigenvalues  yield  the  mutual  information  of  a 
multiple  antenna  radio  channel,  for  which  a  central  limit  theorem  is  proven;  4)  ensembles 
which  yield  robust  estimation  of  a  population  covariance  matrix  with  application  to  array 
signal  processing;  and  5)  a  sample  covariance  matrix  for  which  a  CLT  is  studied  on  linear 
statistics  of  its  eigenvalues,  whose  limiting  empirical  distribution  of  its  eigenvalues  is  stud¬ 
ied  with  application  toward  computing  the  power  of  a  likelihood  ratio  test  for  determining 
the  presence  of  spike  eigenvalues  in  the  population  covariance  matrix.  Details  are  given 
below. 

Information-plus-noise  model.  This  matrix  is  dehned  as 

a  =  {l/N){R^  +  aXr,){Rn  +  aXr,)*, 


where,  R^  and  X^  are  both  nx  N,  the  entries  of  X^  are  independent  standardized  random 
variables,  and  a  >  0.  Here,  both  n  and  N  are  considered  large  and  on  the  same  order 
of  magnitude.  As  mentioned  in  the  previous  interim  reports  the  eigenvalues  of  Cn  are 
signihcant  to  both  array  signal  processing,  where  it  models  the  sample  correlation  matrix 
of  N  “snapshots”  or  recordings  taken  from  a  bank  of  n  antennas  of  signals  (the  columns 
of  Rn)  in  a  noise-hlled  environment  (columns  of  and  in  wireless  communications, 

modeling  schemes  where  the  fading  channel,  {l/\/N){Rn  -|-  crW^),  has  non-zero  means. 
Assume  Hn,  the  empirical  distribution  function  of  the  eigenvalues  of  {l/N)RnR'^,  {Hn{x)  = 
(l/n)(number  of  eigenvalues  of  (1/A^)i?^i?*  <  x))  converges  in  distribution  to  H  a,s  n  ^ 
oo,  and  lim^^oo  n/N  =  c  >  0.  Then  it  has  been  shown  in  reference  [6]  of  the  proposal  that 
Fn,  the  empirical  distribution  function  of  the  eigenvalues  of  Cn,  converges  in  distribution, 
with  probability  one,  to  a  non-random  F,  whose  Stieltjes  transform 

=  f  T~ — dF{X)  >  0 

J  X  —  z 


is  the  unique  solution  for  Az  >  0  to 


(*) 


m  = 


t 


(1  -|-  a‘^cm)z  -|-  cr^(l  —  c) 


dH{t). 


Under  additional  assumptions  on  the  entries  of  Xn  the  following  has  been  shown  in  [1]. 
Assume  [ai,a2]  C  (61,62)  where  61  >  0  and  (61,62)  is  outside  the  support  of  F.  Let 
Cn  =  n/N  and  denote  the  distribution  having  Stieltjes  transform  which  satishes 

(*)  with  c,H  replaced  by  Cn,Hn  respectively.  A  natural  assumption  is  to  impose  (61,62) 
is  also  outside  the  support  of  for  all  n  large.  Then,  it  is  proven  that 


P(no  eigenvalue  of  Cn  appears  in  [ai,  61]  for  all  n  large)  =  1. 


1 


Time  was  devoted  on  showing  the  second  part  of  exact  separation,  namely  the  correct 
nnmber  of  eigenvalnes  on  either  side  of  [ai,6i].  For  Cn  this  shonld  correspond  to  the 
nnmber  of  eigenvalnes  of  {1/ N)RnR^  on  either  side  of  some  interval  (di,  ^2)-  The  strategy 
is  similar  to  that  of  reference  [2]  of  the  proposal,  which  establishes  the  exact  separation 
of  the  eigenvalnes  of  a  class  of  large  dimensional  sample  covariance  matrices.  The  idea 
is  to  systematically  increase  the  nnmber  of  colnmns  of  {Rn  +  uX^),  nntil  the  matrix 
C'^  =  {1/N')[R'^  +  aX'^)[R'^  +  {X'  being  the  nnmber  of  colnmns  of  R'^  +  (^X'^)  is 

close  to  {l/N)RnRn*  +cr‘^I^  I  being  the  nxn  identity  matrix.  This  can  be  achieved  from 
known  resnlts  on  the  extreme  eigenvalnes  of  {1/N')X'^X'^*  and  on  resnlts  established  by 
the  principal  investigator  on  the  cross  terms  {l/N'){R'j^X'^*  +X'^R'^*).  The  increase  in  the 
nnmber  of  colnmns  is  done  in  a  hnite  nnmber  of  steps,  where  at  each  step  no  eigenvalne 
of  C'^  will  cross  the  interval  [ai,  61]. 

The  machinery  in  performing  the  steps  is  all  in  place,  except  for  one  important  factor. 
It  is  crncial  that  there  is  an  interval  in  [ai,bi]  which  does  not  shrink  to  a  point  as  the 
nnmber  of  colnmns  increase,  and  the  principle  investigator  has  not  yet  shown  this  to  be 
trne.  In  reference  [7]  of  the  proposal,  for  each  c^,  Hn,  a  fnnction  Xn{b)  for  6  >  0  is  created 
which  is  associated  with  [ai,  61]  and  increases  on  an  interval,  with  minimnm  and  maximnm 
valnes  being  the  endpoints  of  the  interval  which  contains  [ai,  61]  and  is  ontside  the  snpport 
of  From  simnlations  this  interval  appears  to  be  increasing  in  length  as  c„  decreases 

(corresponding  to  the  addition  of  more  colnmns),  bnt  the  intervals  are  not  simply  nested, 
that  is,  not  necessarily  one  ontside  the  other.  It  is  when  the  principal  investigator  shows 
the  intervals  remain  larger  than  some  positive  length  will  he  be  able  to  complete  the  proof 
of  exact  separation.  When  he  is  able  to  do  so,  he  will  acknowledge  fnnding  from  this  grant. 

Multi-source  power  estimation.  Research  was  conducted  on  a  way  to  estimate 
the  transmit  powers  of  several  sources  in  mult-antenna  fading  channels,  when  the  number 
of  sensors  and  the  number  of  samples  are  sufficiently  large  in  comparison  to  the  number  of 
sources.  Consider  K  stations  which  transmit  data,  where,  for  /c  =  1,...  ,iF,  transmitter 
k  has  Uk  antennas  with  (unknown)  transmission  power  P^.  These  stations  transmit  to  a 
receiver,  a  collection  of  N  sensing  devices.  Let  N  x  Uk  denote  the  channel  matrix 
between  transmitter  k  and  the  receiver.  It  is  typically  assumed  that  the  entries  of  Hk  are 
independent  with  mean  zero  and  variance  (1/A^).  At  time  m  (m  =  1, . . .  ,  M)  transmitter 
k  transmits  vector  G  to  the  receiver,  impaired  by  additive  noise  G 

(cr  >  0).  It  is  also  typically  assumed  that  the  entries  of  and  are  independent 

over  position  and  time,  and  standardized.  At  time  m  the  receiver  senses 

V^kHkxi^^  + 

k=i 


Let  Y  (N  X  M),  Xk  {uk  x  M),  and  W  {N  x  M)  denote  the  respective  matrices  resulting 
from  placing  each  of  as  columns  of  a  matrix.  With  H  =  [Hi, . . .  ,  Hk] 

and  X  =  [Xj" , . . .  ,  X'^]'^ ,  Y  can  be  written  as 

Y  =  HP^Xx  +  aW, 


2 


where,  with  n  =  ni  +  ■  ■  ■  +  n^,  is  n  x  n  diagonal  with  its  hrst  ni  entries  ■,  next 
77-2  entries  ■,  and  so  on. 

The  goal  is  to  estimate  the  powers  Pk  from  the  matrix  Y.  This  has  been  achieved  by 
extending  the  work  of  X.  Mestre  ([2]),  who  has  determined  an  effective  way  of  estimating 
the  popnlation  eigenvalnes  from  those  of  the  sampled  ones  of  the  matrix  described 
in  the  proposal.  This  work  in  tnrn  relied  on  references  [1],  [2],  and  [17]  of  the  proposal. 
Notice  Y  is  the  hrst  N  rows  of 


fHpy^  iN\fx\ 

V  0i  O2  J\wJ^ 


(/tv  N  X  N  identity  matrix,  Oi,  n  x  n,  O2  n  x  N  zero  matrices),  so  that  resnlts  in  [1],  [2], 
wonld  apply,  if  not  for  the  fact  that  the  entries  of  X  and  W  are  typically  not  of  the  same 
distribntion.  In  [3]  the  extension  of  [1]  and  [2]  to  non  iid  entries  is  made,  and  the  following 
is  proven: 

Theorem.  Let  =  {1/M)YY* .  Fix  iF,  Assnme  M  >  N ,  n  <  N ,  M/N  — )■  c, 
each  N/nk  — )■  Cfc,  as  — )■  cx)  and  certain  assnmptions  on  the  size  of  c,  and  the  Cfc’s. 
Let  Xi  denote  the  i-th  smallest  eigenvalne  of  and  s  =  (v)Xi, . . .  ,  VXn)* ■  Then  with 
probability  1  Pfc  — ^  Pfc  as  — )■  cx)  where 


Pk 


NM 

nk{M-N) 


iVi  -  ^7), 

ieMk 


where  A4  =  {A^- ■  ■  ■  ,^-E*=fc+i  TT-i},  the  ?7i’s  are  the  ordered  eigenvalnes  of 
diag(Ai, . . . ,  Atv)  —  (1/A^)'S'S*,  and  the  HiS  are  the  ordered  eigenvalnes  of  diag(Ai, . . . ,  Xn)  — 
{1/M)ss*. 

The  paper  demonstrates,  throngh  simnlation,  the  snperiority  of  this  scheme  over  an¬ 
other  method  of  estimation,  which  nses  moment  methods. 

Mutual  information.  This  work  concentrated  on  establishing  a  central  limit  theo¬ 
rem  on  a  qnantity  fnndamental  to  MIMO  (mnltiple  inpnt-mnltiple  ontpnt)  wireless  com- 
mnnications  systems.  Consider  a  commnnications  system  with  n  transmiting  antennas 
sending  information  to  N  receiving  antennas.  At  an  instant  of  time  data  Xn  G  is 
transmitted  across  the  system,  resnlting  on  the  received  antennas  the  vector  i/nP-^  given 
by 

Vn  Hn^n  T  ^n- 

Here,  i/„,  A^  x  n,  is  the  fading  channel  matrix,  and  Zn  G  is  white  noise  on  the  receiving 
antennas,  having  variance  p  for  each  of  its  components.  In  a  Rician  channel,  it  is  assnmed 
that 

where  X^  is  N  xn  containing  iid  complex  Ganssian  entries,  Rn,  N  x  N ,  T„,  n  x  n,  and  A^, 
N  X  n  are  deterministic  with  Rn  and  Hermitian  nonnegative  dehnite  {rI/‘^  and 
being  their  nonnegative  dehnite  sqnare  roots).  Rn  and  Tn  are,  respectively  the  correlations 
between  the  receive  antennas  and  the  transmit  antennas.  The  An  matrix  allows  for  “direct 
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line  of  sight”  between  each  pair  of  transmit  and  receive  antennas.  According  to  the  classical 
work  of  Shannon,  the  mutual  iuformatiou  of  the  system  is  given  by 

N 

I  =  ^log(p“^Aj  +  1)  =  logdet{p~^Hr,H*  +  /), 
i=i 

where  Ai . . .  ,Ajv  are  the  eigenvalnes  of  I  is  the  rate  at  which,  for  any  small 

probability,  p,  a  coding  scheme  can  be  fonnd  for  which  the  data  can  be  transmitted  with 
decoding  error  eqnal  to  p.  Dne  to  the  entries  of  being  complex  Ganssian,  when  dealing 
with  /  it  is  snfficient  to  assnme  Rn  and  are  both  diagonal.  Consider 

IM  =  i^/N)  logdet(if^if;  +  pi). 

This  qnantity  is  the  mntnal  information  per  receiver  antenna.  It  follows  that 

I  =  NIM -N  log  p. 

Work  on  the  closeness  of  Rnip)  to  a  bonnded  deterministic  valne,  W(p),  determined  by 
Rn,  Tn,  and  as  n  and  N  approach  inhnity,  has  been  done  previonsly  (Hachem,  et.al. 
Ann.  Appl.  Probab.  17(3),  2007,  875-930).  It  is  assnmed  that 

N  N 

(**)  0  <  liminf  —  <  limsnp  —  <  oo. 

n  n  n 

The  flnctnations  of  Tn{p)  are  important  to  nnderstand,  for  example,  in  determining 
the  probability  of  an  outage,  which  is  the  inability  to  transmit  at  a  given  rate.  The  limiting 
distribntional  behavior  of  In{p)  has  been  achieved  in  [4].  In  the  paper  it  is  shown,  when 
the  entries  of  are  iid  standardized  with  hnite  16-th  moment,  and  when  (**)  holds,  there 
exists  a  0^^,  dehned  by  Rn,  T^,  and  A^  ((2.3)  of  the  paper)  for  which 

-X(IJp)  -  EI„(p))  A  N(0, 1), 

V 

as  n  — )■  cx).  The  convergence  of  the  bias  term  N{EXn{p)  —  Vn{p))  is  also  shown  to  converge 
nnder  certain  assnmptions.  In  particnlar,  when  the  entries  of  Xn  are  complex  Ganssian, 
it  is  shown  that 

7V(EX„(p)-I4(p))=0(iV-i). 

This  central  limit  theorem  is  the  hrst  assnming  a  Rician  system,  where  A^  is  non-zero. 
Robust  M-estimators.  Research  was  also  condncted  on  bringing  together  work  on 
large  dimensional  sample  covariance  matrices  and  that  of  robnst  estimation  of  popnlation 
covariance  matrices.  A  large  body  of  work  on  the  latter  has  been  developed  for  the  pnrpose 
of  handling  ontliers  and  heavy-tailed  distribntions.  Let  x  G  be  a  mean  zero  random 
vector  with  covariance  matrice  Cn  —  E{xx*).  The  goal  is  to  acqnire  an  estimate  of  Cn  from 
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N  >  n  independent  samples,  xi, . . .  ,xn  oi  x.  A  robust  M-estimator,  Cn  of  Cn  is  defined 
to  be  a  solntion  to 


* 

'i 


where  -u  >  0  is  a  snitably  chosen  fnnction.  This  matrix,  known  to  exist,  has  displayed 
important  robnst  properties,  and  converges  with  probability  1  to  (7^  as  — )■  oo.  The 
principal  investigator  worked  on  showing  properties  of  Cn  when  n  and  N  both  increase  to 
infinity,  bnt  are  on  the  same  order  of  magnitnde.  This  was  achieved  in  two  cases.  The  first 
is  where  x  =  Xn  =  AnVm  where  An  is  n  x  m  withA^A^*  =  Cn  and  i/n  €  C™'  consists  of 
independent  standardized  entries.  It  was  done  by  comparing  Cn  to  the  sample  covariance 
matrix 


S 


n 


i=l 


The  following  resnlts  are  containing  in  [5].  Assnme  u  :  M"*"  — )■  M"*"  is  nonincreasing,  contin- 
nons,  with  the  fnnction  (^{s)  —  su{s)  nondecreasing,  bonnded  with  snp^  (f){s)  >  1.  Assnme 
the  existence  of  positive  rj  and  a  snch  that  for  all  n  maxj  Edynjl^^"^)  <  a,  and  with 
Cn  =  n/N,  Cn  —  m/n  assnme 


0  <  liminfcn  <  limsnpCn  <  1  and  limsnpcn  <  cxd. 

^  n  n 


Finally,  assnme  all  eigenvalnes  of  Cn  are  contained  in  a  fixed  bonnded  interal  away  from  0 
for  all  n.  Then 

Theorem  1:  (I)  There  exists  a  nniqne  solntion,  Cn  to  (**)  for  all  large  n  a.s.  which 
can  be  obtained  by 

Cn  =  lim 

t— >-oo 


where  =  /  and 


N 


N  \  n 

i=l 


—x*(Z^^^)  ^Xi]xiX 


(II)  With  II  ■  II  denoting  spectral  norm  on  matrices 


U-\l)Cn-L 


0. 


Using  known  properties  on  Sn  Theorem  1  was  applied  to  the  direction  of  arrivals  problem 
in  array  signal  processing.  Assnme  K  signals  are  impinging  on  a  collection  of  n  sensors 
with  angles  9t,  ■  ■  ■  ,9k-  At  each  instant  of  time,  t,  assnme  the  vector  of  data,  xt  G  C”  at 
the  sensors  is  of  the  form 

K 

=  X]  VPks{9k)zk,t  +  crwt, 

k=i 
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where  s(0)  G  is  the  known  nonrandom  nnit  norm  steering  vector  for  signals  impinging 
the  sensors  at  angle  9,  Zk^t  G  C  is  the  signal  sonrce  which  is  standardized  with  hnite  8  +  ry 
absolnte  moment,  iid  across  t  and  independent  across  /c,  pfc  >  0  is  the  transmit  power  of 
sonrce  k,  and  awt  G  is  the  noise  received  at  time  t,  where  w  containing  iid  entries,  stan¬ 
dardized,  with  hnite  8-|-?7  absolnte  moments.  Then  xt  =  AnUt,  with  —  [5'(0)P^/^  cr/], 
5'(0)  =  [s(6'i),...  ,s(6'k)],  P  =  diag(pi,...  ,pk),  and  yt  =  ■  ■  ■  ,ZK,t,wJ)^  G  . 

Assnme  the  above  conditions  on  n,  N,  and  M  =  N  +  K.  It  is  straightforward  to  verify 
that  Cn  has  smallest  eigenvalne  cr^  with  mnltiplicity  n  —  K.  Then  the  following  is  also 
proven. 

Theorem  2:  Let  E  G  be  a  matrix  with  colnmns  containing  orhogonal  eigen¬ 

vectors  of  Cn  corresponding  to  eigenvalne  a^.  Let  denote  the  eigenvector  of  Cn  with 
eigenvalne  Afc,  Ai  <  ■  ■  ■  <  A^.  Then,  as  n  — )■  cx) 


7(^)  -  m  ^  0, 


where 


y(6/)  =  s{eyEE*s{e) 

n 


l=\ 


and 


1  +  ELn- 


En — K 

fc=l 


K+1 

A 


A 


Ai  A/^ 
_  Afc 


j£k 


Ai  Afc 


.^i  Afc 


i  <  n  —  K 
i  >  n  —  K 


with  fii  <■■■<  fin  the  eigenvalnes  of  diag(A)  —  ,  A  =  (Ai,  ■  ■  ■  ,  E)"*"- 

In  the  classic  MUSIC  (mnltiple  signal  classihcation)  algorithm,  if  Cn  were  known,  7 
wonld  be  nsed  to  determine  the  angles  9j,  being  the  zeros  of  7.  With  Theorem  2  a  robnst 
estimate  of  7  in  the  large  dimensional  setting  can  be  made  where  the  estimates  of  the 
angles  are  made  as  the  valnes  which  minimize  7. 

In  [6]  the  above  resnlts  were  extended  to  the  Xi  being  elliptically  distribnted.  Specif¬ 
ically,  it  is  assnmed  that  Xi  —  y/riAnyn,  where  An  is  as  above  except  n  <  m,  the 
are  real  random  with  some  additional  assnmptions,  and  the  yi  are  independent  nnitar- 
ily  invariant  zero-mean  vectors  with  ||yi||  =  m,  independent  of  the  r^.  Assnme  also  that 
(snpg  (/)(s))cn  <  1.  Then  Theorem  1  of  the  previons  paper  holds  trne  in  this  case,  where 
now  (II)  is  now 

||a-*Sn||  ^0, 

where 

1  ^ 

Sn  =  j;^'^VniTan)XiXy 
i=l 

with  ■jn  is  the  nniqne  positive  solntion  to  the  eqnation 


y’njrj'y) 

iV  1  +  CnilniriT) 
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and  Vn{x)  =  Ur,{g^^{x)),  tl^nix)  =  XVr,{x),  gn{x)  =  x/{l  -  Cn(p{x)). 

CLT  of  linear  spectral  statistics  of  a  spiked  population  model.  Consider  the 
sample  covariance  matrix  with  a  spiked  popnlation  model  discnssed  in  the  proposal:  = 

{l/N)Tn^'^XnX*Tn^‘^  where  Tn^^,  n  x  n,  is  the  Hermitian  nonnegative  dehnite  sqnare  root 
of  the  nonnegative  dehnite  matrix  T^,  X^  is  nx  N  consisting  of  i.i.d.  standardized  entries, 
Cn  =  n/A^— )-c>0,  asn— )-cx),  and  for  a  hx  M  the  eigenvalnes  of  are  1  with  mnltiplicity 
n  —  M,  the  remaining  M  eigenvalnes  being  different  than  1.  The  limiting  behavior  of  the 
M  ’’spiked”  eigenvalnes  is  given  in  (2)  of  the  proposal.  Let  Hn  denote  the  empirical 
distribntion  function  of  the  eigenvalnes  of  T„,  and  the  nonrandom  distribntion 

fnnction  corresponding  to  the  empirical  distribntion  fnnction  of  the  eigenvalnes  of  B^,  as 
described  in  the  proposal:  is  dehned  throngh  (x)  =  1—Cn+CnF^^’^^  {x  >  0), 

which  corresponds  to  the  empirical  distribntion  of  the  eigenvalnes  of  {l/N)X*Tri,Xn,  the 
Stieltjes  transform  having  a  nniqne  inverse  given  by 

1  ft 

z{m)  == - hCn  /  ———dHr,{t). 

m  j  1  +  tm 


In  [7]  the  distribntional  behavior  of  linear  statistics  of  eigenvalnes  of  B^  for  arbitrary 
Tn  is  stndied  where  Hn  converges  in  distribntion  to  H.  Let  Ai  >  ■  ■  ■  >  denote  the 
eigenvalnes  of  B„.  When  the  second  and  fonrth  moments  of  Xu  match  that  of  a  real  or 
complex  normal  random  variable,  it  is  shown  that  for  analytic  fnnction  / 


(*  *  *) 


Xn{f]=n(lj2XX)-  j 


f{x)dF 


Cn  yH-n 


converges  in  distribntion  to  a  Ganssian  random  variable  with  known  mean  and  variance 
depending  on  c  and  H.  The  convergence  will  be  the  same  whether  Tn  is  the  identity 
matrix  or  of  spiked  behavior  (M  eigenvalnes  different  than  1),  since  the  limiting  H  will 
not  change.  However,  the  difference  between  the  two  will  be  of  order  1/n  and  in 

(***)  there  is  an  n  mnltiplying  the  difference  between  the  linear  statistic  and  the  centering 
term,  rendering  a  signihcant  difference  between  the  spiked  and  nnspiked  models. 

In  [8]  the  difference  between  the  two  is  derived  np  to  an  0(l/n^)  term.  It  is 

given  as  Theorem  1  in  the  paper.  It  is  applied  to  determine  the  power  of  a  test  introdnced 
in  [9].  Under  Ganssian  assnmptions  and  when  c  <  1  A  corrected  likelihood  ratio  statistic 

L*  =trBn-  log  |Bn  -  n 


denoting  determinant)  is  introdnced  to  test  the  hypothesis 


Hq:  T  =  I  vs.  Hi  :  T^F 
They  prove  in  the  paper  that  nnder  Hq 

L*  -n(^-  ^  log(l  -  Cn)^  N{m{g),v{g)), 


7 


where 

lofff  1  —  c) 

m{g)  = - and  v{g)  =  — 21og(l  —  c)  —  2c. 

2 

This  will  yield  the  probability  of  a  type  I  error. 

Because  of  Theorem  1  in  [8],  the  power  of  the  test  can  be  determined  under  the 
specihc  alternative  hypothesis  Hi'.  has  M  eigenvalues  different  than  1,  having  distinct 
eigenvalues  ai  >  ■  ■  ■  >  ak  with  respective  multipliciites  ni, . . .  ,  n^.  It  is  shown  in  [8]  that, 
under  the  centering  term  f  (x  —  logo:  —  is 


1  ^ 

1  +  -  logtti) 

n  “ 


i=l 


n 


Therefore,  from  (*  *  **),  under 


L*  —n  I  {x  —  log  a:  —  N{m{g),v{g)). 


With  $  denoting  the  standard  distribution  function,  using  signihcance  level  a  (typically 
.05)  it  follows  that  the  asymptotic  power  function  of  the  test  is 

ni{ai  -  1  -  log  g^)  ^ 

A/-21og(l  -  c)  -  2c  ) 


1  -  $  $  ^(1  -  a) 
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